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Abstract. We present the results of a spectroscopic campaign revealing that the non-thermal radio emitter CygOB2#8A is an 
06 + 05.5 binary system. We propose the very first orbital solution indicating a period of about 21.9 days. The system appears 
to be eccentric (0.24 ± 0.04) and is likely seen under a rather low inclination angle. The mass ratio of the components is close 
to unity. The impact of the binarity of this star in the framework of our understanding of non-thermal radio emission from 
early-type stars is briefly discussed. 
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1. Introduction 

Cyg OB2 is one of the most interesting OB association of our 
Galaxy. It is known to harbour a huge number of OB stars, 
among which about 100 O-stars (Knodlseder .2000 . Comeron 
et al. l2002> . Considering its mass, density and size, Knodlseder 
J2000> suggested to re-classify it as a young globular clus- 
ter, the first object of this class in our Galaxy. However, the 
extremely large and non-uniform reddening towards Cyg OB2 
renders its study at optical wavelengths extremely difficult. 
One of the optically brightest O-stars in CygOB2 is 
CygOB2#8A (BD -h40° 4227, my =9.06). This star is classi- 
fied as 05.51(f) (Massey & ThompsonTWi. CygOB2#8A is 
one of the first early-type stars discovered to be an X-ray emit- 
ter with Einstein (Harnden et al. ll979> . Its X-ray emission was 
subsequently investigated with ROSAT (Waldron et al. ll998> 
and A5CA (Kitamoto & Mukai .1996, De Becker|200D- 
A peculiarity of this star is that it is a non-thermal radio emitter 
(Bieging et al. '1989V This non-thermal emission requires the 
existence of a population of relativistic electrons, believed to 
be accelerated in shocks through the Fermi mechanism (Chen 
& White [T994t . An important question is what is the nature of 
the shocks responsible for this acceleration process: intrinsic 
shocks due to radiative instabilities within the wind of a single 
star (e.g. Chen & White 1994), or a wind-wind collision within 
binary systems (e.g. Eichler & Usov J393 1? Concerning the 
latter scenario we note that while the majority of non-thermal 
radio emitting Wolf-Rayet stars are indeed binaries (Dougherty 
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& Williams I2000t . the situation is less clear for O-stars (De 
Becker et al.'2004l). In particular, up to now no evidence for bi- 
narity was put forward for Cyg OB2 #8 A. For instance, Lortet 
et al. J1987> did not detect any bright (Am < 3 mag) visual com- 
panion with an angular separation in the range 0.04 to 1 .5 arc- 
sec. 

In this paper, we present the first spectroscopic monitoring of 
Cyg OB2 #8A revealing that it is a binary system. The obser- 
vations are described in Sect.|2 The spectrum is discussed in 
Sect.|3 and in Sect.l^we present our period determination, as 
well as the first orbital solution proposed for this system. The 
conclusions and prospects are finally given in Sect.|5l 

2. Observations and data reduction 

Spectroscopic observations of CygOB2#8A were collected 
at the Observatoire de Haute-Provence (OHP, France) dur- 
ing four observing runs in September 2000, September 2001, 
September 2002 and October 2003. All observations were 
carried out with the Aurelie spectrograph fed by the 1.52 m 
telescope (Gillet et al. l 19941 1. Aurelie was equipped with a 
2048x1024 CCD EEV 42-20#3, with a pixel size of 13.5 //m 
squared. All spectra were taken with a 600 1/mm grating with 
a reciprocal dispersion of 16 A mm"', allowing us to achieve a 
spectral resolution of about 8000 in the blue range. 

A total of 35 spectra were obtained over the wavelength 
range between about 4455 and 4900 A . The data were reduced 
using the midas software developed at ESO, with the same pro- 
cedure as described by Rauw & De Becker (2004i. Tabled 
gives the journal of our observations, including the radial ve- 
locity measured on the Hei A 4471 line (see Sect. l4.1> . The 
mean signal-to-noise ratio of individual spectra of our data set, 
estimated over a region devoid of lines, is about 250. 
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Table 1. Journal of our observations of Cyg OB2 #8A. The first 
column provides the heliocentric Julian date (in the format HJD 
- 2450000). The second column gives the phase according to 
our orbital solution (see Table|2}. The third and fourth columns 
yield the radial velocities measured on the Hei A 4471 line 
for the two components of the binary system. The last column 
gives the weight assigned to the various RV values in the orbital 
solution. The RVs with a weight superior to 0.5 are those which 
were determined through the simultaneous fit of two Gaussians 
with the simplex method. 



Date 


Phase 


RVi 


RV2 


Weight 






(kms"') 


(kms"') 




1810.466 


0.152 


234 


-74.6 


0.1 


1811.468 


0.198 


55.2 


-97.0 


0.1 


1812.475 


0.244 


36.3 


-110.6 


0.3 


1813.523 


0.291 


54.7 


-84.7 


0.3 


1814.510 


0.336 


59.9 


-111.1 


0.3 


1815.503 


0.382 


56.3 


-121.5 


1.0 


1819.462 


0.563 


34.7 


-86.1 


0.3 


1821.488 


0.655 


10.1 


-58.9 


0.1 


2163.394 


0.261 


394 


-97.2 


0.5 


2164.382 


0.307 


53.1 


-104.4 


1.0 


2165.361 


0.351 


72.3 


-80.5 


1.0 


2167.367 


0.443 


53.7 


-84.4 


0.7 


2168.376 


0.489 


38.7 


-44.5 


0.3 


2169.371 


0.535 


17.0 


-32.6 


0.1 


2170.351 


0.579 


0.7 


-42.9 


0.1 


2170.446 


0.584 


0.5 


-33.0 


0.1 


2518.409 


0.466 


29.9 


-90.7 


0.1 


2520.369 


0.556 


214 


-76.4 


0.1 


2522.361 


0.647 


-1.2 


-44.8 


0.1 


2524.354 


0.738 


-60.8 


34.2 


0.1 


2529.394 


0.968 


-110.5 


65.2 


1.0 


2532.340 


0.103 


6.9 


-54.8 


0.1 


2532.421 


0.106 


94 


-25.5 


0.1 


2533.318 


0.148 


18.8 


-87.1 


0.1 


2533.432 


0.153 


21.1 


-81.5 


0.1 


2919.402 


0.771 


-66.0 


41.2 


0.1 


2922.388 


0.907 


-85.9 


114.4 


1.0 


2923.336 


0.950 


-113.6 


84.2 


1.0 


2925.407 


0.045 


-50.4 


17.9 


0.05 


2926.671 


0.103 


-43.9 


11.1 


0.1 


2928.349 


0.179 


37.0 


-84.4 


0.7 


2929.283 


0.221 


59.0 


-83.1 


0.3 


2930.334 


0.270 


36.6 


-91.4 


0.3 


2934.286 


0.450 


62.2 


-102.0 


0.7 



3. Spectral analysis 

The inspection of our time series of spectra reveals that some 
lines display variations pointing towards Cyg OB2 #8 A being 
a binary system. This is most obvious for the He i A 4471 and 
He II A 4686 lines. The former shows indeed a (partial) deblend- 
ing of lines compatible with a binary system, and the latter un- 
dergoes strong variability even if no clear deblending is ob- 
served. 

We measured the equivalent widths (EWs) of the Hei A 
4471 and Hen /I 4542 lines to determine the spectral types 




Fig. 1. Results of the fit of two Gaussians to the He: A 4471 
(left) and He 11 A 4542 (right) lines for the spectra obtained on 
HJD 2452 164.382 and HJD 2452 922.388 respectively for the 
upper and lower figures. These profiles correspond respectively 
to phases 0.307 and 0.907 ouf our orbital solution. The individ- 
ual components are overplotted. The vertical dashed lines give 
the boundaries of the domain where the fit was performed. The 
lower panels provide the residuals of the fit. 

of both components. To separate the lines of the two stars, 
we performed the measurements on the spectra displaying the 
most pronounced separation between the components (HJD 
2452922.388). We simultaneously fitted two Gaussians to the 
profile, allowing the central position, the width, and the nor- 
malization factor of both Unes to vary. The best fit parame- 
ters were obtained using an N-dimensional downhill simplex 
method (Nelder & Mead 1965 1, with N being the number of 
free parameters (6 in our case). The integration of the two indi- 
vidual Gaussians gives the EWs of the lines of the components 
of the binary. Figure^shows the results obtained for the He i A 
447 1 and He n /1 4542 lines at the date near quadrature selected 
for the spectral type determination. Using the classification cri- 
terion of Mathys ( 1988 1, the EW ratios obtained from these fits 
allow us to infer 06 and 05.5 spectral types respectively for 
the primary and the secondary. We note that spectra obtained 
at the opposite quadrature display similar intensities for both 
components. For instance, the profiles displayed in Fig.[2give 
intensities (as well as widths and consequently EWs) which do 
not diff'er by more than about 5-6%. Consequently, we do not 
find evidence for a significant 'Struve-Sahade' efifect (see e.g. 
Bagnuolo et al. ll999 1. The intensities of the lines in the spectra 
of the two components roughly yield a visual brightness ratio 
of 2 between the primary and the secondary. Since the spectral 
types of both components are similar, this suggests that their 
luminosity classes might be different. With an 06I + 05.5III 
binary system, the large My mentioned by Herrero et al. (I2002> 
is readily explained. 
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4. Cyg OB2#8A as a binary system 

4.1. Period determination 

As a first step, we searched for the period of the system through 
a study of the Une profile variability of our time series. A Time 
Variance Spectrum (TVS, Fullerton et al. 1996 1 computed be- 
tween 4455 and 4890 A on all our spectra reveals a significant 
variability for He i .1 4471, Hen AA 4542,4686, H/3 and Niii 
AA 4634,4641. In each case, the TVS displays a double peaked 
profile suggesting a simultaneous variability of the wings of 
these lines, compatible with a binary scenario. We performed a 
Fourier analysis on our times series using the method described 
by Heck et al. ((1985 1, and used e.g. by Rauw & De Becker 
J2004^ for BD +60° 497. All periodograms display peaks at fre- 
quencies of 0.043 and 0.094 d except for the Hen A 4686 
line which does not show the 0.094 d"' peak. Prewhitening 
of the periodogram (see e.g. Rauw et al.'200H reveals that 
these frequencies are harmonics, corresponding respectively to 
P = 23.26 d and about P/2. The former period is more com- 
patible with our data even though P/2 gives often the strongest 
peak. This could be explained by the typical duration of our 
observing runs (max. 16 d), shorter than the actual period. 

As a second step, we measured the radial velocities of both 
components of the system' using the fitting method described 
in the previous section when the separation of the lines was 
sufficient. When the lines were too heavily blended, we set 
the line widths and intensities to their values determined from 
the fits near quadrature and we used an iterative scheme where 
the line positions were first estimated by eye and subsequently 
changed to improve the residuals. These measurements were 
performed only for the He i A 447 1 line, as the other lines were 
not sufficiently separated. The typical error on the RVs is about 
10 kms"' for the cases where we were able to fit Gaussians^, 
and we estimate that it can reach about 25-35 km s"' at phases 
where the lines are severely blended. The Heck et al. ( I1985> 
Fourier technique was then applied to our RV time series of 
the primary (RVi), the secondary (RV2), and to the difference 
RV1-RV2. The periodograms of the RVi and RV2 series present 
both highest peaks corresponding to periods of 23.3 1 and 21 .88 
d (+ 0.04 d). These values emerge also from the RV1-RV2 data 
set. However, we have folded our RVs with these two periods 
with some success and both values were used as first guesses 
in our search for an orbital solution. 

4.2. Orbital solution 

We obtained the first orbital solution for this system using the 
method described by Sana et al. (2003 1, and used e.g. by Rauw 
& De Becker (I2004t . We assigned various weights to our data 
to take into account the different errors affecting our RV mea- 
surements. Due to the intensity ratio of the lines of the two 
components (~ 2), there was no ambiguity on the identifica- 

' One of our spectra was dropped because of a strong contamination 
of the He i A 447 1 line profile by a cosmic. 

^ This error on the radial velocity corresponds to the standard devi- 
ation determined for the radial velocity of a Diffuse Interstellar Band 
(DIB) at about 4762 A obtained for spectra of similar S/N ratio. 
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Fig. 2. Radial velocity curve of CygOB2#8A for an orbital 
period of 21.908 d. The hexagons (resp. triangles) stand for 
the primary (resp. secondary) RVs. The solid and dashed lines 
yield our best fit orbital solution respectively for the primary 
and the secondary. 

Table 2. Orbital solution for CygOB2#8A. To refers to the 
time of periastron passage, y, K, and a sin / denote respec- 
tively the systemic velocity, the amplitude of the radial velocity 
curve, and the projected separation between the centre of the 
star and the centre of mass of the binary system. Rrl stands for 
the radius of a sphere with a volume equal to that of the Roche 
lobe computed according to the formula of Eggleton ( I1983t . 



Primary Secondary 


P (days) 


21.908 (fixed) 


e 


0.24 ± 0.04 


To (HID-2450000) 


1807.139 ±0.894 


7 (kms"') 


-8.1 ±3.3 -25.0 ±3.6 


A" (km s-') 


82.8 ± 3.5 95.8 ± 4.0 


a sin i (Rg) 


34.8 ± 1.5 40.3 ±1.7 


q = nil/ 1712 


1.16 ± 0.06 


m sin' I (Mo) 


6.4 ± 0.6 5.5 ± 0.5 


Rrl sin i (Rq) 


13.6 ±0.2 14.8 ±0.2 



tion of the lines respectively of the primary and the secondary. 
Table|5] yields the main parameters of the system according to 
our best orbital solution-'. As can be seen from Fig.|5J the sys- 
tem is eccentric with e = 0.24 + 0.04. However, we emphasize 
that the error we obtain on the eccentricity underestimates the 
actual eiTor as a result of the rather heterogeneous phase cov- 
erage of our time series (see Fig.|5}- 

The rather low values of the RV amplitude {K), along with 
the modest values of m sin"' / as compared to the typical masses 
of stars of that spectral type, suggest that the system is seen un- 
der a rather low inclination angle, which would explain why 
the lines remain strongly blended over large parts of the or- 
bital cycle. Consequently, CygOB2#8A is unlikely to display 
eclipses. The projected Roche lobe radii (Rrl sin /, computed 
for an actual orbital separation equal to the semi-major axis) do 
not suggest any mass transfer through Roche lobe overflow. 

' The orbital solution with the period of about 23.3 d turned out to 
be incompatible with a spectrum kindly provided by Dr. A. Herrero 
(private communication). 
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We note that the accuracy of this orbital solution is affected 
by several factors. First, we measured RVs only for the Hei 
A 4471 line, because of the strong blending of all other lines. 
Second, our data do not allow us to state that a simple two- 
Gaussian model is adequate at all orbital phases. Indeed, in- 
dividual lines can present intrinsic variabihty independent of 
their motion due to binarity. In that case, individual line pro- 
files should not be gaussian anymore. Moreover, our data do 
not allow us to reject completely the possibility of a third com- 
ponent that could account for some discrepancies between the 
data and our orbital solution. More data obtained with a larger 
telescope, and with a better spectral resolution are necessary to 
address these issues. 

5. Conclusions and future work 

We present for the first time spectroscopic data revealing that 
CygOB2#8A is an 06 + 05.5 binary system with a period of 
2 1 .9 d, likely seen under a rather low inclination angle. The two 
components are probably of different luminosity classes. The 
system appears to be eccentric, with a mass ratio of the two 
components close to 1 . The binarity of Cyg OB2 #8 A implies 
that the fundamental parameters derived through application of 
a single-star model atmosphere code (Herrero et al. 2002 1 will 
have to be revised. 

The He ii A 4686 line presents a strong line profile variabil- 
ity which will be investigated in a subsequent paper to search 
for the signature of a putative wind-wind interaction. 
The binarity of this system is crucial in the framework of the 
study of non-thermal radio emitters. In fact, 4 out of 8 O-stars 
displaying non-thermal radio emission listed by Rauw r2004l 
are now confirmed binaries, while 2 more are suspected bina- 
ries. This lends further support to the idea that multiplicity and 
hence colliding winds are fundamental ingredients for the ac- 
celeration of the relativistic electrons (Eichler & Usov 1993 1 
that generate the synchrotron radio emission (White Q.985 1. 
ASCA data reveal that the X-ray spectrum of CygOB2#8A 
shows a hard emission component with a temperature of more 
than 15 X 10^ K (De Becker EOHD- Provided that the thermal 
nature of this hard component is confirmed by our forthcom- 
ing observations with XMM-Newton, this temperature could 
be compatible with a plasma heated by a wind-wind collision 
(Stevens et alFiW^. 

However, the short period of this system raises questions 
regarding the detectability of synchrotron radiation in the ra- 
dio domain. The radio photospheres, even at 3.6 cm, are huge 
(more than 1000 and 700 R© respectively for the primary and 
the secondary) and the collision zone of the two winds is deeply 
embedded therein. A similar situation occurs however in the 
case of WR 104 (P ~ 1 yr), where synchrotron radiation origi- 
nates well within the dense wind of the Wolf-Rayet component 
(Monnier et al.l2002V Considering this result, with a mass loss 
rate lowered by a factor of about 5-10 as expected for O stars 
compared to WR stars, even with a shorter separation, we ex- 
pect that a non-thermal radio spectrum could be detected in the 
case of Cyg OB2 #8 A. Information on the filling factors of the 
two winds (if they are clumpy) and on their geometry (if they 
are not spherical) are needed to understand this unexpected 



transparency. Additional studies are foreseen to address this 
issue, including an intense radio monitoring of Cyg OB2 #8 A 
to search for a putative variabihty correlated with our orbital 
period. 

Finally, as a binary system harbouring a small population of 
relativistic electrons, CygOB2#8A may also contribute to the 
yet unidentified EGRET gamma-ray source 3EG J2033-I-4118 
(Benaglia et al. 2001l through inverse Compton scattering of 
photospheric UV photons. This hypothesis will be investigated 
through forthcoming observations of the Cyg OB2 region with 
the INTEGRAL satelhte. 
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